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A vector-dominance two-photon exchange model is proposed to explain the re- 
cently observed production of p°p° and p°(p pairs in e + e~ annihilation at 10.58 GeV 
with the BaBar detector. All the observed features of the data — angular and decay 
distributions, rates — are in agreement with the model. Predictions are made for 
' yet-unobserved final states. 
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£ 1 Introduction 



So far all observed hadronic processes in e + e~ annihilation are in agreement with one- 
photon exchange leading to C = —1 final states. However, the BaBar Collaboration 
has recently presented the first measurement of e + e~ annihilation into pairs of neutral 
vector mesons at 10.58 GeV centre-of-mass energy pQ. The final states p° p° and p° (f) are 
observed with cross sections 

a(e + e-^pp) = (20.7 ± 0.7 stat ± 2.7 syst ) fb, (1) 
a(e + e-^ P 4>) = (5.7 ± 0.5 stat ± 0.8 syst ) fb (2) 

in the centre-of-mass angular range | cos#| < 0.8. These processes characterized by C = 
+1 final states cannot originate from one-photon exchange followed by quark- ant iquark 
fragmentation. 
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An explanation for these processes is that they arise from e + e~ annihilation to two 
virtual photons, with each virtual photon converting to a vector meson. The rates pre- 
dicted by this mechanism can be computed simply and unambiguously using the effective 
vector meson-photon couplings determined from the meson leptonic widths. In this note, 
we present the results for the cross sections and show that they are of the correct size to 
explain the cross section values obtained by BaBar. 



2 The model 

We consider the generic process e + e~ — ► V1V2, where V\ and V2 are neutral C — —I 
vector mesons. It is assumed to proceed with an intermediate e + e~ ^77 process, the 
two photons converting into V\ and V2, with effective couplings ejj\ and e// 2 - 

The cross section in the narrow-width approximation can be simply written in terms 
of the Mandelstam variables s = (k + k') 2 , t = (k — pi) 2 , and u — (k — P2) 2 , where k, k', 
pi, and P2 are the incoming electron, positron, and outgoing vector meson 4- momenta, 
respectively. In the e + e~ centre-of-mass, one gets 



d(T ViV2 _ ( e V ( e V da^ji , 2 
dcos0~U J fj dcos9 ^'"^ 



ml, ml) (3) 



with 



rf(T 7r7 2 * / 2 2x Tia 2 2\^\2{m 2 r +m 2 v )sut + {t 2 + u 2 )(ut-m 2 v m 2 r ) 

, a (my , my ) = j= 1 2 - ^ 1 2 —, 4 

dcosO 1 2 s \/s u 2 t 2 



and where 6 1 is the angle between the electron and the V± momenta, and ~p* the meson 
momentum. The last factor reduces to the familiar t/u + u/t = 2(1 + cos 2 9)/ sin 2 9 
when V\ and V2 are made massless. For my x = my 2 = Mz, Eq. (JHJ) is consistent with the 
well-known expression for the vector part of the e + e~ — > Z Z cross section [2j. 

For masses which are small compared to y/s, the angular distribution differs little from 
the massless case. To leading order in rriy/s, for the VV case, one obtains: 

daw a + cos 2 9 

dcos9 ~ b - cos 2 9 ' ( ^ 

with 

a = 1 + ?^, (6) 

* = 1+4- p) 

For m v = m p and s = (10.58) 2 GeV 2 , a = 1.043 and b = 1.00012. 

For the VV final states involving 2 identical particles, the angular distribution must 
be integrated over only one hemisphere, while for V1V2 (V\ 7^ V2), the integration should 
be done over the full angular distribution. The measurements with BaBar are made from 
in the cos# range from —0.8 to 0.8. 
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To take into account the fact that the p meson cannot be properly described in the 
narrow-width approximation and also to generalize to any vector hadronic final state, it 
is advantageous to rewrite Eq. (J3j) in a more general form: 

d(T VlV2 _ ( « V f dm i f dm l p, 2n d(T HY 2 ,1 2, 

where R{m 2 ) = a(e + e~ — > hadrons)/cr p t. Eq. (jHJ) is a completely general expression for 
the 2-virtual-photon cross section, valid as long as the interference between the products 
of the 2 photons can be neglected. All one has to do is plug in the contribution to R 
for any final state and integrate. The narrow-width limit is obtained by taking R(m 2 ) = 
9-rr/a 2 Y y ee my 5{m 2 — my). 



3 Results 

For the narrow vector mesons the leptonic widths can be directly used to derive the e/fy 
values 

while for the p meson we integrate over the mass distribution, taking as input for R{m 2 ) 
the fit to the annihilation data [HUE]. 

The angular distribution is given in Fig. It is strongly peaked along the beams, 
like the e + e" — > 77 cross section. The behaviour predicted by Eqs. (I3I4I5D is in good 
agreement with the BaBar data pQ. 

The model also allows one to compute the cross sections for mesons with helicity 
A = ±1 and A = 0. As expected the former is dominant, i.e. the vector mesons are 
photon-like with transverse polarisation, leading to a sin 2 9 nt K angular distribution in the 
p — > 7r + 7r~ and — > K + K~ decays. This behaviour is also clearly observed in the BaBar 
data p. 

Integrating numerically Eq. IHlover cos# from —0.8 to +0.8 and over the p mass dis- 
tribution, one obtains 



< p = (21.4 ± 0.7) fb, (10) 



p p 



a™ ^ = (6.15 ± 0.22) fb. (11) 



The BaBar measurements are given with a cut on the p lineshape, retaining only 
events with a mass between 0.5 and 1.1 GeV. The window is likewise defined between 
1.008 and 1.035 GeV. The cross sections in the BaBar conditions are: 



(18.7 ± 0.6) fb, (12) 
(5.3±0.2)fb. (13) 

The agreement with the BaBar results is good, thus confirming the vector-dominance 
two-photon process as the dominant dynamics for V^V® final states. 



th, mass cuts 
°PP 

th, mass cuts 
a p<t> 
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Figure 1: The angular distribution for the processes e + e — ► p p and e + e — > p 4> at 10.58 
GeV as predicted by the vector dominance model in the narrow width approximation. 

The calculation can be straightforwardly applied to other final states, the correspond- 
ing predictions being summarized in Tabled We emphasize that any interference between 
the hadronic products originating from the two virtual photons has been neglected — an 
approximation which is expected to be less valid as the hadronic masses increases on the 
scale of \fs. 

The value obtained for e + e~ — ► J/ip J/ip differs somewhat from the estimate given 
by Bodwin et al [6 . That calculation uses the 2-photon exchange model, with splitting 
of each photon in c c. The authors allow both pairings of the charm and anticharm 
quarks and account for the interference between these diagrams using a non-relativistic 
QCD analysis. They find a cross section of (6.6 ± 3.0) fb, integrating the full angular 
distribution, to be compared with (2.38 ± 0.15) fb for the subset of diagrams that we 
consider. While the new interference terms could be important for J/ip production, since 
2mm is close to y/s, they are not significant for the production of pairs of low- mass states 
such as p and (f). 

The predictions in Table Q beyond the already measured p p and p final states can 
be tested in the future as some of them are within reach of the BaBar and Belle detectors 
with present and foreseen luminosities. 
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Table 1: Cross sections for the production of neutral vector meson pairs V\ V2 in e + e~ 
annihilation at 10.58 GeV. For the final states containing p mesons the results are given 
both with the narrow-width (NW) approximation and the full integration of the mass 
spectrum. The angular integration is restricted to the cos#Vi range from —0.8 to +0.8 
for the results in the last column. The values for the leptonic widths are taken from 
Refs. |S1 IU Ej- The integral over the p mass distribution according to Eq. IHluses a fit to 
the annihilation data [313]. 



v l v 2 


T V2 ^ ee (keV) 


a NW (fb) 


a (fb) 


C"|cos6»|<0.8 (fb) 


p p 


i .U 1 ± U.ii 


izo.oi 


1 1 r ar\ 

Ho. by 


01 on 


OJ 


0.60 ±0.02 


21.11 


20.27 


3.76 


P 


1.27 ±0.04 


33.70 


32.39 


6.15 




5.40 ±0.17 


40.66 


39.61 


10.16 


P *I>(2S) 


2.10±0.12 


13.21 


12.89 


2.60 


P T(LS0 


1.31 ±0.03 


7.83 


7.87 


5.45 


00 




2.23 




0.44 


OJ 




2.83 




0.54 


OJ OJ 




0.89 




0.16 


W 




5.06 




1.46 






2.38 




1.13 


J/ip ip(2S) 




1.49 




0.78 



4 Contributions to R 

The computed cross sections can be extrapolated to lower energies in order to estimate 
the contribution of these C — ±1 two-photon processes in comparison to the dominant 
C = —1 one-photon annihilation into hadrons. The latter one is the needed input to the 
calculations of hadronic vacuum polarization at lowest-order. 

Figure El shows the y/s dependence of the cross sections for the dominant processes, 
computed in the narrow-with approximation and expressed relatively to the pointlike 
cross section (R). The angular integration is still performed from cos#vi — —0.8 to 
±0.8, values which correspond to a typical experimental acceptance. It is found that the 
contribution of these 2-photon exchange processes to the measured R is negligible at all 
energies. In particular, their effect in hadronic vacuum polarization calculations (effect 
which is presently neglected), such as for the anomalous muon magnetic moment or the 
running of a is two orders of magnitude smaller than the present experimental accuracy 
from the R input values. 

Interference effects between C = —1 and C = ±1 amplitudes could occur at the 10~ 2 — 
10~ 3 level, but they cancel for charge-symmetric event detection. So their contribution is 
negligible. 
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Figure 2: The contributions to R from the dominant processes e + e~ — > V\ Vn as predicted 
by the vector dominance model in the narrow-width approximation and integrated in the 
[—0.8 — 0.8] cos 9 range. 

5 Conclusions 

The proposed two-photon exchange vector-dominance model for the processes e + e _ — > 
V\ V2, where V\ and V% are C = — 1 neutral vector bosons, agrees with all the features 
of the BaBar data pQ for the p p and p <fi final states and provides a good description of 
their rates. 

The cross sections for other possible final states, p u, p J/ip, p ip(2S), p T(15), 4> (ft, 
to (j), u to, 4> J/ip, J/ijj J/ip, and J/ip ip(2S) have been predicted. Most of them are within 
reach at the B factories with foreseen luminosities. 

The contamination of these two-photon exchange processes to the measured R ratio 
is much smaller than any foreseeable experimental uncertainty on R. 

Helpful discussions with S. Brodsky and many BaBar colleagues are warmly acknowl- 
edged. 
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